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Abstract. This letter presents generalization of the Landau hydrodynamic solution for 
multiparticle production applied for non-central relativistic heavy ion collisions. Obtained 
solution enables to calculate anisotropy parameter V2 for different rapidities with respect 
to the initial conditions for hydrodynamics. 

PACS: 24.10.Nz - Hydrodynamic models, 25.75.Ag - Global features in relativistic heavy 
ion collisions, 25.75.Ld - Collective flow 

Landau's approximate hydrodynamic solution for particle production in relativistic heavy 
ion reactions pQ [2] is still valid for RHIC energies and reproduce particle spectra for differ- 
ent pseudorapidities well [3]. Experimentally observed azimuthal asymmetries of particle 
production in non-central heavy ion collisions are of high interest nowadays, as it does 
provide more information about the early dynamics of the high-energy nuclear reactions. 
Moreover, RHIC data [7] of elliptic flow for different pseudorapidities shows universal be- 
havior for different nuclei and for different beam energies, which is obtained from the 
model. 

The purpose of this letter is to introduce an analytic solution for relativistic hydrody- 
namic equations, including azymuthally asymmetric transverse expansion. The presented 
solution is approximate, but in comparison to the computational hydrodynamics is ana- 
lytic and transparent. The main assumptions of the model coinsides with original Landau 
approximations and are: i) longitudinal and transverse parts of hydrodynamic equations 
are solved separately; ii) equation of state of ideal relativistic gas, p = e/3, is used to solve 
transport equations; iii) transverse expansion does not include initial flow and is pressure 
gradient driven. 

In this model we solve the equations of energy-momentum conservation: 

(1) = 0, 
where the energy- momentum tensor reads as: 

(2) T» v = (e + p)u»u u -pg» u . 

The solution of local conservation laws (eq. [T]) with the equation of state represents the 
dynamics of the system by relating bulk properties of the matter, such as: energy density, 
e, pressure, p and the four- velocity of the fluid, = u°(l,v). 
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The equations of the hydrodynamic longitudinal expansion in 1+1 dimension, along z 
axis reads as: 



1ZZ 



Q T 00 q T 0z dT 0z g T z 

< 3 > ~w + ih =0 - -w + — = ' 

Solution of the above equations starts by transforming relativistic velocity field components 
to rapidity terms, as: 

(4) u° = cosh rj, u z = sinhr/. 

From the above transformation, naturally follows, that: v z = tanhr/ and (u z ) 2 — (u ) 2 = — 1. 
The details on how to solve the equations of hydrodynamics ([3]) can be found in pQ, [2], [4], 
[5], [6] so here we present only the resit, as we will need it later. For abbreviation following 
variables are used: 

(5) ?? + = l/21n((t + z)/A), 

(6) ) 7 _ = l/21n((t-js)/A) J 

while A is the initial thickness of the system in the beam direction, z. Also, A is the initial 
condition after which equation of state assumed to be valid and evolution equations ([3]) 
are applied. 

The final solution for energy density, e(r/ + ,^_), and rapidity, r](r] + , ??_), is: 

(7) e(r ?+ ,r ? _) = e exp[-4/3(r?+ + ??_ - V r /+ r /-)]> 

(8) rj(ri + ,r]-) = (rj+ -rj-)/2, 

while z = ttanhr/. The above solution of ([3]) is exact and analytic of 1+1-dimensional 
relativistic hydrodynamics, which will be used after solving transverse expansion, to obtain 
the multiplicities of produced particles for different rapidities. 

In order to solve the transverse part of hydrodynamic equations we will follow original 
Landau assumptions with some modifications. For simplicity polar coordinates will be used, 
where four-flow in polar coordinates is: U{ = dxi/dt, uq = (1 — (r 2 + r 2 <fi 2 ))~ 1 / 2 , u r = uqv t 
and energy-momentum tensor ([2]) components are as follows: 

T rr = (e + p) (u°) 2 v 2 + p, T** = {e + p) (n°) 2 t; 2 + p/r 2 

T 0r = (e+p)(u°) 2 v r , = ( e + p)(«°)% 

Assuming that system does not rotate, <j) = = 0, hydrodynamic equation (fTJ) for the 
transverse dynamics at the fixed transverse angle (p is: 



dT 0r dT r 

Inserting energy-momentum tensor expressions to the above equation, we get: 
(10) 4e(n°) 2 ^: + 4 e (n0) 2 ^ + | = 

Following original Landau derivation, the fist term in the above equation is acceleration and 
assumed to be proportional to dv r /dt oc r(t)/t 2 , the second term is set to zero, because v r 
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being comparably small. To simplify the third term Landau used de/dr ~ —e/R, because 
the energy density at the center has value e and is zero at the edge of the system, r = R. 
In the case of peripheral collisions, we do not expect centrally symmetric energy density 
distribution, thus the assumption is modified, as: 

de e(r = R(<j>)) — e(r = 0) 
( ' fr = RlJ) ' 

where R(4>) is a transverse radius of the system, which changes with the angle, as the system 
is not centrally symmetric. We do not know the value of e(R(4>)), so we introduce a new 
function, as f(R(c/>)) = e(r = R(<f>))/e(r = 0), which is a fraction of energy density at the 
edge of the system with respect to the energy density at the center. The function, f(R((/>)), 
naturally must be less than unity for any 4>, as the energy density at the center is higher, 
than at the edge. Now from eq. (fTOj) we express transverse displacement dependance on 
time, as: 

(12) r(t) _(l-/(W 



A{vP) 2 R{(j)) ' 

The last stage of the model is so called conic-flight stage, where energy and entropy fluxes 
stop changing for the fixed cone element litrdr. With the help of the formula (|12|) we 
obtain a hypersurface in space-time, after which hydrodynamics stops and matter streams 
freely towards detectors. Here again we follow original Landau model, assuming a fixed 
transverse distance, r(t) = a, for the conic-flight to start. Thus, we obtain the value for 
the time, when the conic-flight starts, which reads as: 



(13) t\ = 2 cosh m 



aR{4>) 



(1 -/(«))' 



where vP = cosh?y. 

Solution in the conic-flight stage is straightforward, as the energy and entropy does not 
change at a fixed cone element. Connection of transverse and longitudinal solutions is done 
by matching the time t = t\. Knowing, that dS = au°dz at a given time within element dz 
and entropy density, a = ce 3//4 , we express entropy change over rapidity from the energy 
density formula ([7|), as: 

dS 3/4 t 

(14) - = ce exph^ + ^-^T)]^. 

Inserting the solution for time eq. (fT5|) to the entropy equation above and assuming, that 
number of produced particles is directly proportional to the entropy, dN oc dS, one can 
obtain number of particles for different rapidities for a fixed angle 4>. What is still needed, 
is the function f(R((j))), hydrodynamic evolution length, a, in the transverse direction and 
initial starting thickness, A. Besides, one has to normalize the particle distribution (|14j) 
with the total number of produced particles, to obtain the number of produced particles 
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for different rapidities. But, from the definition of elliptic flow, V2- 

J d(f>(dN/d(f>dr]) cos (20) 

(15) V2{7]) = Jd^dN/dW,) ' 

normalization of the total number of particles drops out, yet the dependance on the impact 
parameter b is present. 

In order to show that obtained hydrodynamic solution provides reliable results, we build 
the most simple initial configuration, which is based on Woods-Saxon nuclear density [8] 
and overlap geometry. The only assumption is, that nuclear density is directly propor- 
tional to the energy density. As we do not need the absolute value of energy density, this 
assumption is reasonable enough to obtain qualitative results. 

The thickness function T AA (b,r, 0) of nucleus is obtained, using Wood-Saxon nuclear 
density function [8]: 

n A (r) = 5 — , 

l + exp(^)' 

where R A = 1.12A 1 / 3 - 0.86vl~ 1/3 [fm], d = 0.54 [fm] and n = 0.17/m~ 3 . Integrating 
nuclear density n A (r) over beam axis z one gets it's thickness: 



T l(l ,») = / Jl . < (vS?T?) 



The overlap thickness function with the distance 6 between their centers is T AA (b, x,y) = 
T A (x — 6/2, y) + T A (x + 6/2, y). The transformation to polar coordinates is straightforward 
with x = rcoscj), y = r sin</>. Finally, the function, f(R(<j))), which gives the ratio of energy 
density at the edge of the system with energy density at the center, with fixed 6, reads as: 

(16) f(R(</>)) = T AA (b,R(<f>))/T AA (b,0). 

The radius of the system, R(<f>), is obtained from the geometry on how two circles overlap, 
as: R(<j)) 2 + R((p)b cos <j)+^-R A = 0. 

Merging equations (|13l 1141 [T6|) we connect elliptic flow (|15p with the initial thickness A. 
Now it's possible to use experimentally observed elliptic flow for different pseudorapidities 
[7], which behaves universally for different nuclei collisions (Au-Au and Cu-Cu) and for 
different beam energies (y / s]viv = 64 and 200 GeV). By comparing results with the data, 
we find, that hydrodynamic initial condition, A, must be proportional to the nuclear 
thickness, A = C f (R((j))) / 7 . The following means, that hydrodynamic expansion ([3]) 
starts with azimuthally asymmetric initial thickness, which is wider where initial nuclear 
density is higher. Obtained results are shown in figure ([I]) for Au-Au at 6 = 6fm and 
Cu-Cu at 6 = 3fm reactions for two different beam energies, used at RHIC. We have used 
C = 1.5-Ra and a = R A for all four curves in figure [U 

Quantitative fits to the initial thickness, A, should be done with respect to more realistic 
picture of the energy density profile function, f(R((j))), which might be obtained from 
the distribution of participating nucleons and number of binary collisions with respect 
to the nuclear cross section, at the same time including beam energy dependance to the 
transverse dynamics. These more sophisticated models for the initial configuration are 
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Figure 1. Elliptic flow dependance on rapidity for Au-Au at b = 6fm 
and Cu-Cu at b = 3fm for two different collision energies. Beam rapidity 

Ubeam = '« N N /m N ) . 

in development, but for now we can conclude, that Landau hydrodynamic solution and 
it's assumptions works not only for particle multiplicity spectra, but for elliptic flow and 
elliptic flow scaling as well. 

The generalized solution of Landau hydrodynamics incorporates transverse asymmetries 
into initial configuration and can be compared with wider amount of the experimental 
data. In this way giving deeper insight into the early dynamics of the system. The validity 
criteria of the obtained solution are no different from the original Landau hydrodynamic 
solution. Using just a toy model, we have found that initial thickness of the system is not 
a flat disk, but is proportional to the azimuthally asymmetric nuclear density. 
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